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In supersymmetric theories, cosmic strings produced in the early Universe often 
have a width of TeV scale, while the tension is much larger. In a scaling regime, an 
I infinite cosmic string releases significant fraction of its energy in the form of string 

loops. These thick string loops lose their energies efficiently by particle emissions, 
and hence it may have effects on cosmological observations. We study cosmological 
implications of string loops with TeV scale width in detail and derive constraints 
^ ' on the tension of the string. Implications on future gravitational wave detectors are 

H ' also discussed. 



1 Introduction 



Although the cosmological evolution scenario in the early Universe well before the big- 
bang nucleosynthesis (BBN) is not known yet, it is expected that the early Universe 
experienced some stages of spontaneous symmetry breaking (SSB) as the cosmic temper- 
ature decreases. In association with the SSB, cosmic strings may be formed P depending 
on the topology of the vacuum manifold [2]. Besides, cosmic superstrings emerge at the 
end of the brane inflation [31 H]. Cosmic (super)strings, if detected, will provide a way to 
search for the early Universe and high-energy physics. Their properties and cosmological 
consequences have been widely considered in many papers. 

On the other hand, supersymmetry (SUSY) [5] is one of the promising candidates of 
the physics beyond the standard model, since it provides a natural solution to the gauge 
hierarchy problem. SUSY opens up a window toward the energy scale of the grand unified 
theory (GUT), and hence we can consistently discuss the physics of the early Universe in 
the framework of SUSY. Formations of cosmic strings in the framework of SUSY GUT 
were discussed in the literature [0]. Interestingly, a SSB which produces cosmic strings 
may be associated with SUSY breaking. Suppose that there is a flat direction (p in the 
scalar field space, which is ubiquitous in tree-level supersymmetric theories, and that it 
has the soft SUSY-breaking negative mass-squared — m^|0|^ around the origin. Then, 
the vacua deviate far from the origin of the field space and the scalar fields have large 
nonzero vacuum expectation values (VEVs). The width of the cosmic string is roughly 
determined by the curvature of the scalar potential at the origin and the tension is roughly 
given by the square of the VEV |7j. Since the energy scale of the soft SUSY breaking 
term, m, is of (9(TeV), the width of the cosmic string is estimated as ~ (TeV)~^, unlike 
the ordinary GUT string whose width is expected to be about (lO^^GeV)"^. The VEV, 
on the other hand, is determined by the balance between the negative soft mass term 
and higher dimensional operator, which is suppressed by the cutoff scale M, which can 
be as large as the Planck mass Mp. In this case two mass-dimensional parameters which 
appear in the potential of the flat direction have extremely large hierarchy. The VEV can 
be larger than m by many orders of magnitude. This may push up the tension of cosmic 
strings to a cosmologically relevant value. 

A similar situation has been considered in the context of thermal inflation [SI [U [TOl [H] . 
Thermal inflation, which is the short-lasting inflationary period introduced as a solution 
to the cosmological moduli problem [121 [IS] ? requires the scalar field (called flaton) with 
light mass and large VEV in order to get a sufficient amount of e-folds to dilute away the 
moduli. It is often assumed that the flaton have negative soft mass term at the origin, 
and is stabilized at a large field value due to non-renormalizable terms. If it has a U(l) 
symmetry that is broken due to the flaton VEV, cosmic strings are produced at the end 
of thermal inflation with properties just described above. 

The cosmic strings with TeV-scale width and much higher scale tension have char- 
acteristic cosmological effects beyond which ordinary strings have, such as gravitational 
waves (GWs), anisotropy of cosmic microwave background (CMB) and so on. In particu- 
lar, their extreme thickness enhances the particle production at a cusp on a loop. A cusp 
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is a spiky and extremely Lorentz-boosted region on a cosmic string loop which appears 
several times per one period of the oscillation of the loop. When a cusp appears, two 
parts on the both sides of the cusp overlaps each other and this induces nonperturbative 
production of the scalar particles [HI [15] . This effect clearly becomes more efficient as the 
width of the string get thicker. The produced scalar particles decay to lighter particles 
in the standard model sector. If this happens during the BBN epoch, it can disastrously 
modify the hght element abundances. Besides, the decay products may contain stable 
particles, which may be a candidate of dark matter (DM), and may provide too much 
matter abundance to be consistent with current observations [HI . These cosmological 
effects from the cusp annihilation crucially depends on the typical loop size created per 
Hubble time, parameterized by a, which is highly uncertain up to now. These observa- 
tions motivate us to study cosmological effects of TeV-width strings in detail, and obtain 
constraints on the parameter space of the loop size a and the tension /i through GWs, 
DM and BBN. Cosmic strings associated with a flat direction were considered in several 
works [m dni [201 [21 [22] and especially Ref. [22] partly referred to the cosmological con- 
straint from such cosmic strings. In this paper, we further investigate detailed constraints 
to a and G/i of cosmic strings with TeV scale width from various cosmological effects and 
also argue detectability of GWs emitted from these strings at future experiments. 

This paper is organized as follows. In Sec. [21 a simple SUSY model for the spontaneous 
f/(l) breaking is presented and show that cosmic strings with TeV scale width naturally 
appears. In Sec. [3] evolution of string loops are briefly shown taking account of the GW 
emission and particle emission from cusps. In Sec. [Hwe discuss how the particle emission 
and GWs from string loops affect cosmological observations. In Sec. show constraints 
on the string tension and typical loop size, and also model parameters. We conclude in 
Sec. [a 



2 Cosmic string models with TeV scale width 

Let us consider superfields (p-^ and 0_, which have the charge -|-1 and —1 under the 
additional U{1) gauge symmetry, respectively. One of the candidates of an additional 
f/(l) is the [/(1)b-l symmetry, but we do not specify a concrete setup. We also assume 
that there is R-symmetry and both and 0_ have R-charge +l/n (n is a positive 
integer). These assumptions prohibit any renormalizable terms including in the 

superpotential and the allowed non-renormalizable term is 

w = (1) 

where M is the cut-off scale of the low-energy effective theory. The scalar potential is 
given by 

V = Vf + Vd + V^oft, (2) 



where 

Vr. = 



vf = Ti^i0+i"i0-r"-'(i0+r + i0-r), (3) 



2 



Vo-'^(\4>,?-\4>-\'?, (4) 

Ko,. = -ml\4,,? - ml\^4' - (-4^^ + h.c.) . (5) 

where g is the couphng of the additional U{1) and soft masses m+ and m_ are expected 
to be of order of TeV. A is also of order of TeV and we have implicitly redefined scalar 
fields such that A is real and positive. Here, we use same letters for superfields as their 
scalar components. We assume that the soft mass terms for 0+ and (j)- is negative. The 
D-term potential forces scalar fields to the D-flat direction, = 10-1- The F-term and 
soft SUSY breaking terms produce the global minimum in this direction alQ 



A + \ A^ + (4n - 2){m\ + ml 



l/(2n-2) 



(6) 



An -2 

arg0+ = — arg0_. (7) 
Parametrically, 

where v is the scale of the VEVs of the scalar fields and m is the scale of the coefficients 
of soft SUSY breaking terms, that is, m ~ m+ ~ m„ ~ A. Since the vacuum manifold is 
S^, this potential leads to the emergence of cosmic strings if the U{1) is restored at first 
and then spontaneously broken after that. The width of cosmic string w is roughly given 
by the inverse of the curvature of the scalar potential at the origin as 

w ~ ~ (TeV)"^ (9) 

The tension /i is roughly given by the square of the VEV of the scalar field, namely 

_ (mM2"-3)i/{«-i). (10) 

The feature ^/Ji ^ m does not depend on the detailed model construction. It is rather 
a generic feature as long as the soft SUSY breaking mass triggers the SSB. If we take 
M ~ Mp and m ~ TeV, ~ {W^GeVf for n = 2 and ~ {W^GeYf for n = 3. These 
correspond to G^x ~ 10~^^ and G^x ~ 10~^°, respectively. Although these values are much 
smaller than the current upper limit (G/i < 10~^), cosmic strings with such a tension still 
deserve cosmo logical interest, as discussed below. 

The potential naturally causes a short-period of secondary infiation, called thermal 
infiation [SI HI [TOl [11] . During infiation, and 0_ get masses of order of the Hubble 
parameter through supergravity effect, 

^H = i^fnf(c4-|0+P + C-|0-P), (11) 



^Strictly speaking, the true global minimum shifts from the D-flat direction if ^ to_. However, 
the soft term is subdominant to the D-term around ([7]), then the shift is small. 
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where c+ and c_ are 0{1) constants and assumed to be positive here, and i^inf is the 
Hubble parameter during inflation. If Hi^f > m, which is satisfied in most inflation 
models, 0+ and 0_ are trapped at the origin during inflation. After inflation, the energy 
stored in the inflaton is converted to the thermal bath through the reheating process. At 
the origin, the U{1) is restored and hence the U{1) gauge filed is massless. Then there 
appears finite temperature corrections to the scalar potential as 

W = ctT2(|0+P + |0-P). (12) 

Here ct is a 0{1) positive coefficient. This stabilizes 0+ and 0_ at the origin for T ^m. 
When the temperature falls below Tbe ~ y/mv, the potential energy of the scalar fields 
(~ m^f ^) exceeds the energy of the thermal bath and thermal inflation starts. During this 
period the preexisting radiation and matter is diluted exponentially and its temperature 
goes down. The unwanted relics such as the gravitino and moduli are diluted away, 
and hence it is appealing for solving the cosmological moduli problem [121 113] ■ When 
it falls down to Tend ~ n^', the finite temperature correction becomes negligible and the 
true potential minimum appears. Then the scalar field rolls down toward the potential 
minimum and thermal inflation ends. After thermal inflation, cosmic strings are formed 
since the U{1) symmetry is spontaneously broken by the VEVs of 0±. 

The flaton decays after thermal inflation, and the radiation dominated universe restarts. 
The reheating temperature after thermal inflation must exceed a few MeV for the suc- 
cessful BBN [23] • The decay rate of the flaton is estimated as 

r = 7—, 13 

V 

where 7 is a numerical constant and the reheating temperature is given by 

1/4 



^ lOOMeV (^) f ^) ( ' ) (^) , (14) 

where gt, is the number of relativistic degrees of freedom at T = Trh- For example, the 
flaton couples to the Higgses as [21] 

W = X^H^H,, (15) 

IVl p 

yielding a right magnitude of the higgsino mass of fin = Xv^/Mp. In this case, the flaton 
decays into Higgs boson pair and we obtain 7 ~ (yU///m)^ and hence 7 takes rather wide 
range of values depending on fin- For m = 1 TeV and 7 = 0.1, requiring Trh ^ 0(MeV) 
leads to V < lO^^GeV and we need v < 10^^ GeV for m = 100 GeV and 7 = 0.1. This 
condition sets another upper bound for the tension of cosmic strings. 
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3 Decay of cosmic string loops 



Once formed, the distribution of cosmic strings in the universe obeys the so-called scaling 
law. In the scaling regime, a few infinite cosmic strings exist per Hubble horizon. In order 
for the strings to fall into the scaling regime, cosmic strings must collide and reconnect 
with each other and cast their energies into the string loops. Cosmic string loops lose their 
energies by several processes, which may leave characteristic signatures on cosmological 
observations. We consider two mechanisms in which cosmic string loops lose their energies; 
GW emission and particle emission. 

(i) Gravitational wave emission 

GW emission from cosmic string loops [25] is usually thought of as the main process 
of the energy losses of loops. GW bursts are emitted mainly from cusps and kinks 
(discontinuous inflections on loops) and cusps are more efficient sources compared 
with kinks. The energy emission rate of GWs of frequency uj from a cusp obeys a 
power law, dE/du oc (x!~'^^^, and the total energy loss rate is found by summing up 
the emission rates of all modes as 



where F ~ 50 is a numerical constant. 

(ii) Particle emission 

When a cusp is formed on a loop, two string branches on each side of the cusp 
overlap. It is expected that scalar particles are created due to non-perturbative 
effects. Since it is difficult to estimate precisely the non-perturbative effect, we 
here assume that the conversion from the energy of the condensate into particles is 
efficient enough to approximate the total energy of the created particles to be that 
stored in the overlap region. Then the energy loss rate of a loop through particle 
emission is estimated as [U] (see also Appendix), 



where p is the number of time when cusps appear on a loop in one period of its 
oscillation, w is its width and I is its circumference. Note that particle emission 
becomes more important for thicker or smaller loops. This is simply because the 
overlap region near a cusp becomes relatively larger for such loops. The exponents 
of fi in f|T6|) and (|T7|) imply that for smaller Gfi particle emission becomes more 
important compared with GW emission. 

The loop size at which the efficiencies of these effects are equal is given by 



(16) 




(17) 




(18) 
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For / < Z=, a loop loses its energy mainly by particle emissions. 

Cosmic string loops shrink by loosing their energy and finally disappear. The loop 
length evolves according to 



dl 



This can be solved as 



t-U 



k-l 
1= 




2 A^- 



(19) 



(20) 



where U is the time at the birth of the loop and li is the initial loop size. Here we make 
an assumption that the initial loop size created at the cosmic time ti is given by U = ati, 
with a being a constant smaller than unity. The magnitude of a remains unknown. Its 
suggested value ranges from 0.1 EZ] to some powers of Gfi [2S1 EH]- Intermediate or 
mixture results are shown in Refs. [301 EB 1321 ES] • Considering these uncertainties, we 
will treat a as a free parameter. 

Let us consider some limiting cases. If / ^ /=, GW emission dominates the energy 
loss and then we have I = —TGfi. Therefore, loops become shorter at the constant rate : 



l-li = -TGfi{t - ti) for / > 



(21) 



This is also directly obtained from (|20l) . If a = Zj/tj ^ TGfi, a loop disappears soon after 
it is born and if a ^ TGfi, it survives much more than one Hubble time. On the other 
hand, if Z <^ Z=, we obtain I = — Pa/j- Therefore, we obtain 



3/2 



3/2 



2L 



(t - ti) for / < L 



(22) 



The fractional change of the loop length in one Hubble time at the birth is roughly given 
by 

Pl\l- (23) 



and this becomes unity when 



~1~ 



/V / 

1/3 



(24) 



Therefore, if a ^ (p'^w/ti)^^^ a loop disappears within one Hubble time and if a ^ 
[wp^ /ti)^/'^ it takes longer than one Hubble time for a loop to disappear. 

Particle emissions from cusps extract extra energy from loops, then the lifetime of a 
loop becomes shorter. This alters the distribution of loops in the universe and as a result 
interesting cosmological signatures may be implied as shown below. 
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4 Effects of cosmic string loops on the universe 



4.1 Abundance of cosmic string loops 

As explained, we parametrize the typical length of the string loop / by a constant a@such 
that the loop length is aU at the production time ti. The cosmic string network reaches 
the scaling regime, where there is 0{1) infinite strings in each Hubble horizon. In this 
regime, each infinite string abandons its large portion in the form of loops in each Hubble 
horizon per one Hubble time. Therefore, the number density of loops at t which are born 
between ti and ti + dti is given by 

^"^^'^^U^a-Vf^Vcit., (25) 



dti * \ a{t) 

at the cosmic time t. Here, (a(tj)/a(t))^ represents the dilution due to the cosmic expan- 
sion. The length of a loop at the time t is found by solving (120|) . 

4.2 Gravitational waves 

GW bursts from cusps on loops overlap each other and form the stochastic GW back- 
ground. The spectrum of GW background was calculated in Refs. |25l [SH |35l [361 EZl [38] . 
In the calculation of the GW background, we should be careful that the GW background 
consists of only bursts which overlap other ones. When the contributions of many bursts 
to the background are summed up, we must omit bursts which come to the observer solely 
(which we dub as "rare bursts"). This was pointed out in Refs. [251 [361 [37] ^-^id we follow 
the formalism of Ref . [37] . 

The amplitude of the GW background without omitting rare bursts is given by the 
integral 

arM!)-^ffi^jms.^.i)'^. (26) 

where / is the frequency of the GW, z is the redshift at which the burst occurs, / is the 
loop size, /i(/, /) is the present value of the strain of the GW emitted by the loop with 
size / at redshift and d^R/dzdl x dzdl is the burst rate which reach the observer emitted 
by loops with size / ~ / -|- at redshift z ^ z + dz. h{f, z, I) is given by [25] 

where C ~ 2.68 and r{z) = ao j^^_^~^a{t)^^dt is the proper distance. d'^R/dzdl is derived 
as follows. The rate of bursts emitted at redshift z z + dz hy loops which are born at 
ti ti + dti is 



d^R , , 1^2 2p dndV 



dzdt, = -ei———^—-— ——Q{l{t{z),ti))dzdti, (28) 



dzdti A ""{1 + z)l{t{z),ti)dti dz 
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frequency[Hz] 

Figure 1: The spectrum of the GW background, ^cwif), for various values of G/i with 
a = 0.1. Sohd hues take into account the effect of particle emission and dotted ones do 
not. 



where 9m = ((1 + z)fl)^^^^ is the extent of the burst, dV/dz x dz is the proper spatial 
volume between the redshifts z and z + dz^ l(t,ti) is given as a solution of ( l20l) and 6 
is the step function. We get d'^R/dzdl by converting variables from {z,ti) to (z, /) using 

m 



d^R d'^R 



dzdl dzdh 




TGfi + a 



1 



1 + v^TA 



(29) 



Rare bursts are omitted from the sum in the following way. Since stronger GWs are rarer, 
it is plausible to set the cut-off strain /i^, above which GWs are thought of as rare bursts. 
Then, we set from 

d'^R 



R 



dh 



h. 



dz- 



dhdz 



(30) 



where z^, is the redshift at which cosmic strings appear. Eq. (l30l) means that GWs of larger 
strain than h^, are observed less than once in a period of GWs themselves f~^. GWs of 
smaller strain than h^, are observed many times during the period and overlap each 
other. Therefore, the stochastic GW background spectrum is estimated as 



47r^ 
3^ 



dhh^ 



°° , d'^R 
dz- 







dhdz 



(31) 
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Figure 2: The spectrum of the GW background, flcwif): for various values of a with 
Gfi = 10~^. Sohd hues take into account the effect of particle emission and dotted ones 
do not. 

We show the spectrum of the GW background for various values of Gfi with a = 0.1 
in Fig. [H and for various values of a with Gfi = 10~^ in Fig. [2l Other parameters 
are taken to be w = (lTeV)~^,p = l,r = 50. The solid lines take into account the 
effect of particle emission. Dotted ones do not include it, which correspond to w = 0. 
The shapes of the solid lines are identical to those shown in previous papers [351 l37] . 
GWs of higher frequencies are emitted earlier. The plateau in the high frequency region 
corresponds to GWs emitted in the radiation-dominated era and the downward-sloping 
region corresponds to those from the matter-dominated era. In the high frequency region, 
it is seen that ^gw calculated including particle emission is suppressed compared with 
that neglecting it. This is because such high frequency GWs are emitted at the time 
when the size of loops is so small that the energy of loops is converted into particles 
more efficiently. The suppression is milder than shown in Ref. [22], since in Ref. [22] it 
is assumed that when a loop becomes smaller than /= it instantly disappears converting 
its whole energy to particles. On the other hand, we here trace the evolution of loop size 
even after particle emission dominates GW emission, and such a loop contributes to the 
GW background until it completely disappears. 

We see in Fig. [2] that reducing a shifts the spectrum toward higher frequency direction 
and suppresses the amplitude. The reason why the whole spectrum goes right is simply 
that loops become small and unable to emit low frequency GWs. The lower cut-off of the 
frequency of the GW background from loops is given by ~ (ato)~^ where to is the present 
age of the universe. GWs of frequencies u ~ (ato)~^ are emitted by loops born within 
one present Hubble time. When a is extremely small, there are no GWs of low frequency 
emitted by loops, and the GW background in the low frequency region consists of GWs 
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from kinks on infinite strings [H] . The amplitude in the plateau becomes larger as a gets 
larger by the following reason. Loops emitted from infinite strings behave like matter, 
and their energy density decreases proportional to a^^. Therefore, the ratio of the energy 
of loops born in the radiation-dominated era to that of the radiation increases as time 
goes by until loops completely decay. As a result, the energy fraction of GWs from loops 
become larger. 

4.3 Big-bang nucleosynthesis 

String loops are continuously produced from infinite strings during BBN epoch. It is 
expected that scalar particles emitted from cusps on loops soon decay into lighter par- 
ticles producing many hadrons, since the lifetime of the particle must be shorter than 
1 sec in order for successful reheating, as already explained. Since they emit high energy 
particles, the standard BBN prediction on the light element abundances may significantly 
be affected. Constraints on the energy injection from BBN were studied in Ref. HO] 
and the most stringent bound comes from overproduction of primordial D or ^Li due to 
hadro-dissociation of ''He, which leads to Ap^is/s < lO^^^GeV at T ~ 10 keV, where 
Apvis is the injected energy density per Hubble time and s is the entropy density. 

The energy injection from cusps in one Hubble time around T ~ lOkeV is estimated 
as 

Apvis ^ / dti — -f^Jjr: rT0(^(^BBN, i*)) x tsBN, (32) 

Jo "f^i V H'-BBN, ti) 

where ^bbn is the cosmic time at T = 10 keV and dn{t-BB^,U) / dti is given by Eq. f l25|) . 
Thus BBN provides a constraint on the parameters of string loops. 

Let us focus on some limiting cases. We can classify properties of string loops de- 
pending on whether GW emission dominates particle emission or not, and whether loops 
survive more than one Hubble time or not. The loops disappear as soon as they are 
produced for a <^ TGjj, if GW emission dominates, and for a <^ {p^w/tY^^ if particle 
emission dominates, for the loops created at t. On the other hand, the GW emission is 
dominant for VGjJi > p{w/atY^'^ if loops soon disappear, and for FG/i > p{tp/wY^^ if 
loops are long lived for loops disappearing at t. These are summarized in Fig. [31 

(1) Short-lived and GW emission dominated loops : This corresponds to the upper- 
left region in Fig. |3l The energy density of the string loops is comparable to that of the 
infinite strings with a scaling regime, 

Ploop(t) ^ Gpprit). (33) 

The visible particle emission within a Hubble time per entropy density is estimated as 

^PVmit) ^ Ploop(t) EpE _ PrW / A^A .g^x 

s s s yV^atJ 

(2) Short-lived and particle emission dominated loops : This corresponds to the lower- 
left region in Fig. [31 The energy density of the string loop is given by Eq. (!33|) . Since the 
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all of the loop energy goes into the visible sector, we obtain 



s s 

(3) Long-lived and GW emission dominated loops : This corresponds to the upper- 
right region in Fig. |3l The main contribution to the loop energy density at the time t 
comes from those disappearing at t, which were created at ti ~ TG^t/a. Thus the energy 
density of the string loops exceeds that of the infinite strings with a scaling regime, 

Ploop(t) ^ Ploop(ti) - GuPrit) (^Y^ ■ (36) 



Therefore, the visible particle emission from these loops is estimated as 

Apvis(t) ^ pioop(t) EpEt ^ pr(t) (aw)^/^ 
s ~ s pl{t;ti) ~ s r2G/itV2- 



(37) 



(4) Long-lived and particle emission dominated loops : This corresponds to the lower- 
right region in Fig. [31 The dominant contribution to the loop energy density comes from 
those with the length of l{t;ti) ~ {p'^t'^wY^^, which were born at ti ~ (p'^t'^wY^^ /a. The 
energy density of loops is given by 



Ploop(i) - Ploop(ti) - Gpprit) (^^^ • (38) 

The visible particle emission is then estimated as 

1 /6 

Apvis(t) _ Pioop(t) EpEt _ pioop(t) _ P^^'^Xq^ f (39) 



s s pl{t; ti) s s \p'^w 

From these expression, it is evident that the visible particle emission is efficient at earlier 
epoch. The BBN constraint, however, is most stringent at T ~ 10 keV and hence particle 
emission from loops are constrained from BBN dominantly at around T ~ 10 keV. 
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a 



Figure 3: Parameter regions where loops lose energies mainly through GW emissions 
and particle emissions, and where loops are short-lived and long-lived for w = (lTeV)~^. 
Thick-solid (thin-dashed) lines correspond to loops created at T = IkeV (IGeV). 



4.4 Dark matter production 

It is expected that scalar particles emitted from cusps decay into standard model particles 
and SUSY particles with roughly equal branching ratio as long as the latter process is 
kinematically allowed. At least two lightest supersymmetric particles (LSPs) are produced 
per 0-decay into SUSY particle pair. If the i?-parity is conserved, the LSP is stable and 
contributes to the measured DM abundance. This may explain DM or the condition 
that the density of this decay products must not exceed the current DM abundance sets 
constraints on the loop parameters. 

Denoting the fraction of the energy that is converted into DM particles to the total 
energy emitted from cusps by e, the LSP abundance is given by 

w I ait] \ _ , , , . , 

^LSP = — / dU dte—{t,U fiJ—— e tco - tfo , 40 

Per Ju Jtt, dti Y l{t,ti) \a{to)J 

where is the cosmic time when cosmic strings appear, tfo is the cosmic time when the 
LSP freezes out and pcr is the critical density. LSP freeze-out occurs when the temperature 
drops to Tfo ~ lOGeV. The upper limit of the integral for t, the cosmic time when LSPs 
are emitted, is given by 



iUti) = < ' ' 7 (41) 
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where tco(^j) is defined by l{tco{ti),ti) = 0. In the integral the contribution from LSPs 
emitted before the freeze-out is omitted since they thermahze soon. As described in 
Ref. [T7j, Eq. lHOj) is rewritten as 



1 dn 
^LSP = — / dti2e—{to,ti)iJ,L 

Per Jt, 




hoik) 



eikomeiUu)), (42) 



where 



and 




for ti < tfo 
for ti > tfo 

for tcoiti) < to 

l{to,ti) for tcoiti) > to. 



We expect that 



e ~ 



(43) 



(44) 



(45) 



where mLSP is the mass of a LSP and m,^ is that of the scalar particle of which cos- 
mic strings consist. Hereafter, we set mLs? = lOOGeV, m<^ = ITeV and e = 0.1. 
The condition that ^lsp calculated by P2|) does not exceed the current observed value 
f^DM^^ — 0.1123 IH] leads to the constraint on a and G/i. 

Analytic expressions for ^lsp in some limiting cases are similar to those in the previous 
subsection, except that we need to multiply an additional factor e to them and that they 
should be evaluated at the earliest possible epoch after DM decouples from thermal bath. 
This is because the contribution on the DM abundance goes larger for loops created at 
earlier epochs [17]. Therefore, the requirement that Q-lsp should not exceed the current 
observed value limits the energy injection from cosmic strings around the LSP freeze-out, 
in contrast to the BBN constraint. 

There is a subtlety which deserves to be commented about. If the LSP freezes out in the 
radiation-dominated era, the redshift at the freeze-out is simply given by l + Zfo = Tfo/Tg, 
where Tq is the today's temperature of the universe, and tfo is simply obtained. Analytic 
expressions for I^lsp becomes similar to those in the previous subsection. However, if the 
reheating temperature after the thermal infiation, which is given by f lT^ . is lower than Tfo, 
LSP freeze-out occurs while the oscillation of fiaton dominates the universe. As explained 
in |l2], even during the oscillation of fiaton, there is radiation whose temperature drops 
as T oc a~^/^. Therefore, the redshift and the time at LSP freeze-out is given by 



l + Zfo = 



fo 



8/3 



RH, 



^1 + 2;r,h), ifo ^ 



' 1 + 2:rh \ 



3/2 



tRR, 



(46) 



where 1 + zrh = Trh/^o and tRH are the redshift and the time at the completion of 
reheating, respectively. The dominant contribution on the DM abundance again comes 
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from LSPs emitted around LSP freeze-out. However, whether loops are short-hved or 
long-hved, the energy of LSPs at time t is roughly given by 



ApLSP (t ) ~ Gmot it) , (47) 

where ptot is the total energy density of the universe. This is because the energy density 
of loops is not enhanced even if they are long-lived since both loops and oscillating fiaton 
decrease their energy in proportion to a~^. 

4.5 CMB anisotropy 

It is well known that the cosmic string network induces the CMB anisotropy. However, 
observations such as WMAP have shown that cosmic strings cannot be the main source of 
the CMB anisotropy and hence the string contribution to the cosmic density fluctuation 
is constrained. According to the recent calculation, this bound reads G/i ^ 5 x 10^'' 
This is independent of the parameter a. Hereafter, we concentrate on the case where G/i 
is smaller than 10~^. 



5 Cosmological constraints 

Now let us derive the constraints on the tension (Gp), and loop size (a), and also find the 
parameter region where future experiments are sensitive. We take account of constraints 
from BBN and LSP production as described in the previous section. We include future 
or on-going GW detectors such as ultimate DECIGO BBO, BBO-correlated [IB] , 
LISA [17| and advanced LIGO [IHj as the limits and sensitivities on GW emission. We 
also use the current pulsar timing limit jUJ |50] and future sensitivities from SKA [51] on 
the GW background from cosmic strings. The sensitivity curves of future GW experiments 
are found in Ref. [52j . 

The result is shown in Fig. HI Solid lines take into account the effect of particle 
emission in the calculation of the GW background, while dotted lines do not. The solid 
and dotted lines corresponding to sensitivities of LISA and SKA and the current limit 
from the pulsar timing overlap each other, so dotted ones can be hardly seen. 

First, we comment on the curves corresponding to GW experiments. The shapes of 
the curves reflect the behavior of the spectrum when the values of a or Gp are varied as 
described in Sec. 14.21 Some curves are horizontal for small a. This is because loops cannot 
emit GWs detectable by the corresponding experiments and only the GW background 
produced by kinks on infinite strings contibute. As already mentioned, for smaller Gp, 
the effect of particle emission is more crucial. It may be apparently strange that in 
the curves of BBO and LIGO, larger a region is more sensitive to the effect of particle 
emission, since particle emission is more important for smaller loops. However, for smaller 
a, GWs of frequency to which each detector is sensitive were emitted later from larger 
loops, as shown in Fig. [2l 
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a 



Figure 4: Parameter regions excluded by cosmological constraints and probed by future 
experiments on a-G^x plane. Colored regions are excluded by current observational limits 
or cosmological constraints. The area above each line can be probed by corresponding 
future experiments. Solid lines take into account the effect of particle emission into 
calculation of the GW background, while dotted lines do not. 



The excluded region from BBN is separated into two areas. The triangular region 
appearing in the right side of Fig. H] is in the region where loops are long-lived. In the 
upper half of it GW emission dominates over particle emission (case (3) in subsection 14. 3p . 
while in the lower half the reverse holds true (case (4)). In the quadrangular region at 
the upper-left corner, loops are short-lived. Near the right edge GW emission is dominant 
(case (1)) and near the lower edge particle emission is dominant (case (2)). 

As explained before, the constraint from LSP overproduction is similar to that from 
BBN in that both constrain the visible energy injection. However, the former is more 
sensitive to the energy injection at the LSP freezeout. Since the LSP freezeout takes 
place much earlier than BBN, the loop size is smaller and particle productions tend to be 
dominant at that time. This excludes the large parameter region as shown in Fig. |H On 
the boundary of the prohibited region, the horizontal part is in the region where loops are 
short-lived and particle emission dominated. On the other hand, on the downward- sloping 
part they are long-lived, while particle emission still dominates. Note that in Fig. HJ the 
constraint from LSP overproduction is obtained assuming the LSP freezes out in the 
radiation dominated era. If we use f ll4l) as the reheating temperature, the lower boundary 
of the prohibited region by LSP overproduction becomes nearly horizontal over whole 
range of a, since LSP freeze-out occurs during oscillation of flaton and particle emission 
dominates over GW emission for the parameters near the lower boundary (it can be seen 
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Figure 5: Parameter regions excluded by cosmological constraints and probed by future 
experiments on a-M plane in the case of n = 2. The meaning of the lines are same as 
those in Fig. HI 



from (jUj) that the energy of LSPs emitted at freeze-out depends on only when loops 
are particle emission dominated). It is interesting that the tension is much more severely 
constrained than the bound from CMB (G/i < 5 x 10~^). This limits models of SSB with 
a soft SUSY mass term producing cosmic strings. Furthermore, the BBN and DM bounds 
put considerably severer constraints than the current pulsar timing experiments. 

The constraints in Fig. H] can be translated into those on the parameters M in Eq. ((Tj) 
and a. The results for n = 2 or n = 3 are shown in Fig. [5] and [6l respectively. As discussed 
in Sec. [21 it is difficult to raise the tension as high as the GUT scale keeping the width to 
be TeV scale for n = 2. Therefore, cosmic strings in the model of tt. = 2 can be searched 
only by future space-borne interferometers as long as M is at most the Planck scale. On 
the other hand, for n = 3, the tension of cosmic strings can be as large as the GUT scale, 
and hence various experiments and cosmological bounds set stringent constraints on M 
and a. 

6 Conclusion 

In this paper, we have investigated cosmological constraints on cosmic strings with TeV- 
scale width. We have seen that this type of strings naturally appears when the scalar 
field is the fiat direction and it has a soft SUSY breaking negative mass-squared term, 
since the scale of the soft mass m is (9(TeV). Loops arising from such cosmic strings lose 
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Figure 6: The same as Fig. Hlbut for n = 3. 



the significant part of their energies through particle emission from cusps, in addition to 
GW emission. We have considered effects of the particle production on BBN and DM 
abundance, and the stochastic GW background produced by loops including the effect of 
particle production into the evolution of loops. We have derived constraints on the tension 
fi and the loop size a from BBN, DM, and current pulsar timing limits and also found the 
parameter region which can be probed by future experiments. From the result, we have 
seen that the effects of particle emission on the GW spectrum are important for future 
GW experiments and for the BBN and DM bounds for the cosmic strings with TeV-scale 
width. Specifically, future space-based GW detectors cover very wide parameter ranges 
of interest. 



Appendix 

A Energy loss from a cusp 

Let us write the string trajectory as x((t, t) where cr is the spacelike parameter on the 
string. The most general solution to the Nambu-Goto action is given by [2] 

x(a, t) = ^ {a(or -t)+h{cr + t)}, (48) 

with constraint |a'| = |b'| = 1. The cusp is defined as the point where a' + b' = 0. The 
coordinate is taken so that x = 0, cr = at the cusp. Near the cusp, it can be expanded 
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as 

x(a) ^ ix'^a^ + ^x-a^ (49) 

where the subscript means the value evaluated at a = 0. Its time derivative is written 
as 

i(^) ^ li-^ + b;,) + ^(-a;; + K)a + i(-a- + K)a'. (50) 

From the condition |a'| = |b'| = 1, we soon obtain slq^q = 0, lagp + slq^q = and 
similarly for bo- Using them, we obtain |xp ~ 1 — Ixq^ct^. Thus the the Lorentz factor 
near the cusp is estimated as 



Let us estimate the overlap region near the cusp. Taking into account the Lorentz 
contraction of the string, its radius is given by ~ w/'-f, where w is the string width. 
Notice that the cusp direction is parallel to Xg and orthogonal direction, which coincides 
with the direction of the cusp motion, is parallel to Xq'. Thus the overlap region is given 
by — o"c ^ ^ CTc, where ac = [w/(7|xo |)]^/^. Substituting |xo| ~ l/l and |xo'| ~ 1//^, we 
obtain ac ~ y/vJl, where / is the loop radius. Assuming that the field annihilates efficiently 
and emit particles in the overlap region, the total energy emitted in the form of visible 
particles at the cusp is estimated as 

^cusp = /iCTc ^ Ai(w/)^/^ (52) 
Since the frequency of cusp formation is ~ l/l, the power radiated by the cusp is -Ecusp — 
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